Gauge-Invariant Formalism with Dirac-mode Expansion 
for Confinement and Chiral Symmetry Breaking 
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We develop a manifestly gauge-covariant expansion of the QCD operator such as the Wilson loop, 
using the eigen-mode of the QCD Dirac operator p = 'y'^D'^. With this method, we perform a 
direct analysis of the correlation between confinement and chiral symmetry breaking in lattice QCD 
Monte Carlo calculation on 6* at /3=5.6. As a remarkable fact, the confinement force is almost 
unchanged even after removing the low-lying Dirac modes, which are responsible to chiral symmetry 
breaking. This indicates that one-to-one correspondence does not hold for between confinement and 
chiral symmetry breaking in QCD. In this analysis, we carefully amputate only the "essence of 
chiral symmetry breaking" by cutting off the low-lying Dirac modes, and can artificially realize the 
"confined but chiral restored situation" in QCD. 

PACS numbers: 12.38.Aw, 12.38.Gc, 14.70.Dj 



I. INTRODUCTION 

Nowadays, quantum chromodynamics (QCD) has been 
established as the fundamental gauge theory of the strong 
interaction. However, nonperturbative properties of low- 
energy QCD such as color confinement and chiral symme- 
try breaking [l| are not yet well understood, which gives 
one of the most difficult problems in theoretical physics. 
The nonperturbative QCD has been studied in lattice 
QCD and various analytical frameworks 

In particular, it is rather interesting and important to 
examine the correlation between confinement and chiral 
symmetry breaking p^ - [l8| . since the direct relation is 
not yet shown between them in QCD. The strong cor- 
relation between them has been suggested by the almost 
simultaneous phase transitions of deconfinement and chi- 
ral restoration in lattice QCD both at finite temperature 

[l^ and in a small- volume box [H . 

The close relation between confinement and chiral sym- 
metry breaking has been also suggested in terms of 
the monopole degrees of freedom |12| - [l4j . Here, the 
monopole topologically appears in QCD by taking the 



maximally Abelian (MA) gauge 2^24|. For example, 
by removing the monopoles in the MA gauge, confine- 
ment and chiral symmetry breaking are simultaneously 
lost in lattice QCD [l^, fl3. (The instantons also disap- 
pear without monopoles j22|.) This indicates an impor- 
tant role of the monopole to both confinement and chiral 
symmetry breaking, and these two nonperturbative QCD 
phenomena seem to be related via the monopole. 

However, as a possibility, removing the monopoles may 
be "too fatal" for most nonperturbative properties. If 
this is the case, nonperturbative QCD phenomena are 
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simultaneously lost by their cut. 

In fact, if only the relevant ingredient of chiral sym- 
metry breaking is carefully removed, how will be confine- 
ment? 

In this paper, to get the answer, we perform a direct 
investigation between color confinement and chiral sym- 
metry breaking in lattice QCD, using the Dirac-mode 
expansion in a gauge-invariant manner [25|. 



II. GAUGE-INVARIANT FORMALISM WITH 
DIRAC-MODE EXPANSION 

We newly develop a manifestly gauge-covariant expan- 
sion of the QCD operator such as the Wilson loop, using 
the eigen-mode of the QCD Dirac operator p = 
and investigate the relation between confinement and chi- 
ral symmetry breaking. 



A. Gauge Covariant Expansion in QCD instead of 
Fourier Expansion 

In the previous studies [2^, [l^ , we investigated the rel- 
evant gluon-momentum region for confinement in lattice 
QCD, and found that the string tension tr, i.e., the con- 
fining force, is almost unchanged even after removing the 
high- momentum gluon component above l.SGeV in the 
Landau gauge. In fact, the confinement property origi- 
nates from the low-momentum gluon component below 
l.SGeV, which is the upper limit to contribute to a. 

The previous study on the relevant gluonic modes was 
based on the Fourier expansion, i.e., the eigen-mode ex- 
pansion of the momentum operator p^^. Because of the 
commutable nature of [p^,p'^] = 0, all the momentum 
p^ can be simultaneously diagonalized. which is one of 
the strong merits of the Fourier expansion. Also it keeps 
Lorentz covariance. 
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However, the Fourier expansion does not keep gauge 
invariance in gauge theories. Therefore, for the use of 
the Fourier expansion in QCD, one has to select a suit- 
able gauge such as the Landau gauge [26l - [28j . where 
the gauge-field fluctuation is strongly suppressed in Eu- 
clidean QCD. 

As a next challenge, we consider a gauge-invariant 
method, using a gauge-covariant expansion in QCD in- 
stead of the Fourier expansion. In fact, we consider a 
generalization of the Fourier expansion or an alternative 
expansion with keeping the gauge symmetry. 

A straight generalization is to use the covariant deriva- 
tive operator D^^ instead of the derivative operator 
9^. However, due to the non-commutable nature of 
[D'^,D'^] 7^ 0, one cannot diagonalizc all the covariant 
derivative (/i ~ 1,2,3,4) simultaneously, but only 
one of them can be diagonalized. For example, the eigen- 
mode expansion of D4 keeps gauge covariance and is 
rather interesting, but this type of the expansion in- 
evitably breaks the Lorentz covariance. Then, we con- 
sider the eigen-modc expansion of the Dirac operator 
p — ^f^D^^ or — D'^D'^, since such an expansion 
keeps both gauge symmetry and Lorentz covariance. 

In particular, the Dirac-mode expansion is rather in- 
teresting because the Dirac operator p directly connects 
with chiral symmetry breaking via the Banks-Casher re- 
lation Q and its zero modes are related to the topo- 
logical charge via the Atiyah-Singer index theorem |29| . 
Here, we mainly consider the manifestly gauge-invariant 
new method using the Dirac-mode expansion. Thus, the 
Dirac-mode expansion has some important merits. 

• The Dirac-mode expansion method manifestly 
keeps both gauge and Lorentz invariance. 

• Each QCD phenomenon can be directly investi- 
gated in terms of chiral symmetry breaking. 



In this paper, we adopt the hermite definition of the 
7-matrix, 7^ = Then, p is anti- hermite and satisfies 



y,x 



(2) 



The normalized eigen-state \n) of the Dirac operator p 
is introduced as 



p\n) = iXn\n) 



(3) 



with A„ G R. Because of {75,-^} = 0, the state 75|n) is 
also an cigcn-state of p with the eigenvalue —iXn- The 
Dirac eigenfunction 



(4) 



obeys pipn{x) = iXnipnix), and its explicit form of the 
eigenvalue equation in lattice QCD is 

4 

^ ^lt^P^{x)'lpn{x + fl) - U'_^(.t)V>„(x - A)] 

= iXn'ipnix). (5) 

The Dirac eigenfunction ipnix) can be numerically ob- 
tained in lattice QCD, besides a phase factor. 

According to ?7^(a;) — >■ V{x)Uf_i{x)V\x + fi), the gauge 
transformation of ipnix) is found to be 



Ipnix) ^ Vix)-ipnix), 



(6) 



which is the same as that of the quark field. To be strict, 
for the Dirac eigenfunction, there appears an irrelevant 
ri-depcndent global phase factor e^'^'^^^\ according to the 
arbitrariness of the definition of ipn{x). 

It is notable that the quark condensate (qq) , the order 
parameter of chiral symmetry breaking, is given by the 
zero-eigenvalue density piO) of the Dirac operator, via 
the Banks-Casher relation 



B. Eigen-mode of Dirac Operator in Lattice QCD 



lim lim 7rp(0). 



(7) 



Now, we consider the Dirac operator and its eigen- 
modes in lattice QCD formalism with spacing a in the 
Euclidean metric. On the lattice, each site is labeled by 
X = (xi, a;2, X3, 2:4) with being an integer. In lattice 
QCD, the gauge field is described by the link-variable 
Uf,{x) = e"^9^^(^) 6 SU(iVc), where g is the QCD gauge 
coupling and A^i(a;) G su(iVc) corresponds to the gluon 
field. 

In lattice QCD, the Dirac operator p = j^D/^i is ex- 
pressed with Ufj^ix) as 

^"^^y = o;; XI [C^^(a;)(5x+/i,y - U^f,ix)S^-f,J , (1) 



where the convenient notation U-^ix) = Uj^ix — fi) is 
used. Here, Ci denotes the unit vector on the lattice in 



/i-dircction 



Here, the spectral density of the Dirac operator is defined 

by 



P(A) 



V ^ 



('5(A-A„)), 



(8) 



with the four-dimensional volume V . Also, the zero- 
mode number asymmetry of the Dirac operator p is 
equal to the topological charge (the instanton number) 

Q = y|;^ / d^x Tr (G^i/G^i/), which is known as the 
Atiyah-Singer index theorem, Index(^)=(5 [29| . 



C. Kogut-Susskind Technique to reduce Dirac 
Eigenvalue Calculation 

In calculating the eigenvalue of the Dirac operator p^ 
we use the Kogut-Susskind (KS) technique [3, which 
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is often used to remove the redundant doublers of lattice 
fermions. Here, the use of the KS technique is just the 
practical reason to reduce the calculation of the Dirac 
eigenvalues. In fact, the result of the Dirac-mode projec- 
tion, which will be shown in Sec. IV, is unchanged, when 
the Dirac operator is directly diagonalized. 
In the KS technique, one introduces 

T(x)^7r^72^73-^7r, (9) 

which satisfies T{x — (i) = T{x + /t), due to 7^ = 1 for 
each jjL. As a key relation, "diagonalization" of all the 
gamma matrices 7^ can be done as 

T{x)^^^T{x + fi)^ 7^^{x)l (10) 

with the staggered phase rj^{x) defined by 

^^{x) ^ 1, 7^^{x) = (-l)-i+-+-.- (^ > 2). (11) 

By introducing Xn{x) = T'' {x)ipn{x), the eigenvalue 
equation is reduced into 

4 

^ 'nf^i^)P^,{x)Xn{x + A) - U-f,{x)Xn{x - //)] 

= iKXnix), (12) 

which does not include spinor indices. In the KS method, 
the Dirac operator ^^D^^ is replaced by the KS Dirac op- 
erator rj^Df^ with dropping the spinor degrees of freedom, 
which reduces the lattice- fcrmion species from 16 to 4 Q. 

To summarize the KS method, there is no spin indices, 
and the explicit form of the KS Dirac operator 77^!?^ is 

{v^iD^),^^y = V,^{x)[Uf,{x)6.^+ii,y - U^^{x)6.^-f,^yll3) 

The spinless cigenfunction Xn{x) satisfies Eg. dT^ . i.e., 

4 

VtiDfiXnix) = ^^(77^£)^):r,yX«(2/) = iKx7i{x)- (14) 

y M=i 

In the KS method, the chiral partner 'j54'n{x) reduces 
into (— l)^i+^^+^^+^''Xn(a^), which is an cigenfunction of 
with the eigenvalue —iXn. 
Using the KS technique 0, [11 , the Dirac-mode number 
L'^ X NcX 4 is reduced to be x Nc on the lattice. 
The actual number of the independent Dirac eigenvalue 
A„ is about x Nc/2, due to the chiral property of the 
Dirac operator, i.e., pairwise appearance of ±An, apart 
from chiral zero modes. 



III. OPERATOR FORMALISM IN LATTICE 
QCD 

To keep the gauge symmetry, careful treatments are 
necessary, since naive approximations may break the 



gauge symmetry. Here, we take the "operator formal- 
ism" [23, as explained below. 

We define the link- variable operator U±fj, by the matrix 
element of 

{x\U±i_,\y) = U±i_,{x)Sx±ii,y (15) 

Note that [/^ and [/_^ are hermite conjugate as the op- 
erator in the Hilbert space in the sense that 

{y\Ul\x) = Ul{y)5y+^^, = Ulix - /i)<5..-A,y 

= U-^{x)5x^f,^y ^ {x\U-^\y). (16) 

In the operator formalism, Eq.([5l) for the Dirac eigen- 
state is simply expressed as 



^Y.^^.{U^.-U-^^\n)=^K\n). (17) 



In the KS method, where the spinor degrees of freedom 
is dropped off, one identifies Xn{x) = {x\n), and then 
Ea. ([T^ for the KS Dirac eigen-state is expressed as 



— Yfj^{U,,-U-^,)\n)^iK\n), (18) 



where f}^ is defined by the matrix element of {x\f}^\y) ~ 
Vti{x)5x,y Owing to ri^{x ± fi) = 77p(x), one finds 
VfiU±fi = U±^t)^, so that there is no ordering uncertainty 
in the KS Dirac operator in Eg. (fT8| . 

In the following, we mainly use the ordinary Dirac op- 
erator jfj.Df^ and the spinor cigenfunction ipn{x) = {x\n). 
When the KS method is applied, one only has to use 
the identification of Xn{x) = {x\n) in the following ar- 
guments. The final results are the same between both 
calculations based on jfj.D^ and rj^j^D^. 

The Wilson-loop operator W is defined as the product 
of C/^ along a rectangular loop, 

JV 

W^X{U,,^lJ,fi,,---lJ,^. (19) 

For arbitrary loops, one finds X^aLi ~ ^- "aote 
that the functional trace of the Wilson-loop operator W 
is proportional to the ordinary vacuum expectation value 
(W) of the Wilson loop: 
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tr ^ {X1\U^^\X2){X2\U^.^\X3){X3\U^JX4) ■ ■ ■ {xN\Uf^^\xi) 

Xi,X2,--- ,Xfi 

2 N 

tr^{x\Uf,^\x + fii}{x + fiilU^^lx + ^ fik) ■■■{x + ^\U,,j^\x) 



k=l 



k=l 
N 



tT{Uf,, ix)Uf,., (x + fii)Uf,^ {x + J2M--- U^,^ (x + Y^ Afe)} 

X 

(W) ■ Tr 1. 



k=l 



k=l 



Here, "Tr" denotes the functional trace, and "tr" the 
trace over SU(3) color index. 

The Dirac-mode matrix element of the link-variable 
operator J7^t can be expressed with ipn (x) : 

{■m\U\n) = 'y^^{m\x){x\Uf^\x + jj.){x + fj.\n) 

X 

= Yi,l{x)U^{x)^^{x + fL). (21) 



r 



(20) 



Although the total number of the matrix element is very 
huge, the matrix element is calculable and gauge invari- 
ant, apart from an irrelevant phase factor. Using the 
gauge transformation ^ , we find the gauge transforma- 
tion of the matrix element as 



X 

^ • y{x)U^.{x)V\x + A) • V{x + li)i:n{x + A) 

X 

= Y ^V<^)V iL{x)^n{x + A) = {m\Uy\n). 



(22) 



To be strict, there appears an n-dependcnt global phase 
factor, corresponding to the arbitrariness of the phase 
in the basis |n). However, this phase factor cancels as 
g-itpngJVi. — \ between \n) and and does not appear 
for QCD physical quantities including the Wilson loop. 

In the practical lattice-QCD calculation, we adopt the 
KS technique to reduce the computational complexity, 
as mentioned in Sec H-C. In the KS method, instead 
of '!/'„(x), we use the spinless eigenfunction Xn(x) of 
the KS Dirac operator rjf^Df^, with the identification of 
Xnix) = {x\n)^ and the KS-reduced matrix element of 
is expressed as 

{m\U\n) = 'y^^{m\x){x\U^\x + fl){x + fi\n) 

X 

= Y^li^^U^i^)Xnix + fi). (23) 

X 

In the arguments in the next section, the same results are 
obtained between the calculations based on the original 



Dirac operator 7^7^^ and the KS Dirac operator rj^D^. 

IV. DIRAC-MODE EXPANSION AND 
PROJECTION 

From the completeness of the Dirac-mode basis, 
E„ W){n\ 1, we get 

= EEH("^I^I")("I (24) 

m n 

for arbitrary operators. Based on this relation, the 
Dirac-mode expansion and projection can be defined (25| . 
We define the projection operator P which restricts the 
Dirac-mode space, 

P^^|n)(n|, (25) 

where A denotes arbitrary set of Dirac modes. In P, the 
arbitrary phase cancels between |n) and {n\. One finds 
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= P and = P. The typical projections are IR-cut general, but it would not be important for the argument 
and UV-cut of the Dirac modes: of large-distance properties such as confinement. From 

the Wilson-loop operator 

Firs E ^uv^ E (26) 

lA„l>AiR |A„|<Auv 

Using the projection operator P, we define the Dirac- 
modc projected link-variable operator, 



N 



(28) 



k=l 



ij^ ^ PU^P =J2J2 \m)HU^\n}{r 



(27) 



During this projection, there appears some nonlocality in 
I 



we define the Dirac-mode projected Wilson-loop opera- 
tor, 



JV 



k=l 



(29) 



Then, we obtain the functional trace of the Dirac-mode projected Wilson-loop operator 

N 

Tr W"" = Tr [] U^^ = Tr J/^^^f/;, • • • JJ^^ = Tr PU,,PU^,P- ■ ■ PU^^P 



k=l 



tr E {ni\U^,^n2){n2\Uy,2Wz) ■ ■ ■ {nN\Utj.MWi)^ 



(30) 



r 



which is manifestly gauge invariant. Here, the arbitrary 
phase factor cancels between \nk) and (nfc|. Its gauge 
invariance is also numerically checked in the lattice QCD 
Monte Carlo calculation. 

From Tr W^{R,T) corresponding to the R x T rect- 
angular loop, we define the Dirac-mode projected inter- 
quark potential, 

V^iR) ^ - hm iln{Tr M^^(i?, T)}, (31) 

which is also manifestly gauge- invariant. In the unpro- 
jected case of P = 1, the ordinary inter-quark potential 
is obtained apart from an irrelevant constant, 

V{R) = - lim ;i;ln{Tr M^(P,T)} 



irrelevant const., (32) 



T~>oo T 

- hm iln(W^(i?,r)) 

T— s-oo I 



because of Tr W = (W) ■ Tr 1, as was derived in Eg. ([201 



V. ANALYSIS OF CONFINEMENT IN TERMS 
OF DIRAC MODES IN QCD 

We consider various projection space A in the Dirac- 
mode space, e.g., IR-cut or UV-cut of Dirac modes. With 



this Dirac-mode expansion and projection formalism, we 
calculate the Dirac-mode projected inter-quark potential 
V^{R) in a gauge- invariant manner. In particular, using 
IR-cut of the Dirac modes, we directly investigate the role 
of low- lying Dirac modes to confinement, since the low- 
lying modes are responsible to chiral symmetry breaking. 

As a technical difficulty of this formalism, we have 
to deal with huge dimensional matrices and their prod- 
ucts. Actually, the total matrix dimension of {m\Ufj_\n) 
is (Dirac-mode number)^. On the lattice, the Dirac- 
mode number is x NcX 4, which can be reduced to 
be X Nc, using the KS technique 0, [l], as mentioned 
in Sec.II-C. The actual number of the independent Dirac 
eigenvalue A„ is about L'^ x Nc/2, due to the chiral prop- 
erty of p, i.e., pairwise appearance of ±A„. Even for 
the projected operators, where the Dirac-mode space is 
restricted, the matrix is generally still huge. In addi- 
tion, we have to deal with the product of the huge ma- 
trices (mlU^ln) in calculating the Wilson loop. Then, 
at present, we use a small-size lattice in the numerical 
calculation. 

In this paper, we use SU(3) lattice QCD with the stan- 
dard plaquette action at /3 = 5.6 on 6^ at the quenched 
level. The lattice QCD Monte Carlo calculation is done 
using the pseudo-heat-bath algorithm. The lattice spac- 
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ing a is estimated as a ~ 0.25fm, i.e., 
which leads to the string tension a ^ 
the inter-quark potential. Then, the 



a-i ~ O.SGeV, 
0.89GeV/fm in 
total volume is 



V = (6a)^ ~ (l.Sfm)'*, and the momentum cutoff is 
TT/a ~ 2.5GeV. 

On the 6^ lattice, the Dirac-mode number is G"' x3x4 = 
15, 552, which is reduced to be 6^ x 3 = 3, 888 using the 
KS technique. In fact, the KS Dirac operator ?7^I?^ and 
the KS-reduced matrix element (m|[/^|n) are expressed 
by 3,888 x 3,888 matrix. Considering the pairwise ap- 
pearance of A„ and — A„ apart from chiral zero modes, 
the actual number of the independent Dirac eigenvalue 
A„ is reduce to be around G"' x 3/2 = 1,944. 

To diagonalize the KS Dirac operator ?7^D^, we use 
LAPACK [iOl- For the statistical error on the lattice 
data, we adopt the jackknife error estimate (H. 

We show in Fig. 1(a) the spectral density p(A) of the 
QCD Dirac operator p. The chiral property of p leads 
to p(— A) = p(A). Figure 1(b) is the IR-cut Dirac spectral 
density 



Pir(A) = p{X)e{\X\ - Air) 



with the IR-cutoff Air = 0.5a-i ~ 0.4GeV 
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(33) 
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FIG. 1: (a) The spectral density p(A) of the Dirac oper- 
ator in lattice QCD at /3=5.6 and 6"*. The negative-A re- 
gion is omitted, because of p( — A) = p(A). The volume V is 
multiplied, (b) The IR-cut Dirac spectral density pir(A) = 
p(A)6i(|A| - Air) with the IR-cutoff Air = 0.5a"^ ~ 0.4GeV. 



Note that, using the eigenvalue A„; the quark conden- 



sate {qq) is obtained as 



E 

A„>0 



m 
2m 



A^ -I- m? 



V 



(34) 



where v is the total number of the zero mode of Tf). Here, 
apart from the zero mode, the eigenvalues appear as pair- 
wise, which makes {qq) real. (In lattice QCD, one has to 
take account of the doubler contribution, which can be 
regarded as flavor at the quenched level.) Then, in the 
presence of the IR cut Air for the Dirac eigen-mode, the 
quark condensate is obtained as 



(w)Air 



2m 



(35) 



We show in Fig. 2 the lattice QCD result of the quark 
condensate (99)air as the function of the current quark 
mass m in the presence of IR cut Air . 




FIG. 2: The lattice QCD result of the quark condensate 
{<jq) as the function of the current quark mass m in the 
presence of IR cut Air = 0.5, 1.0, 1.5[a"^]. The vertical axis IS 
normalized by the original value of (qq) without cut. A large 
reduction is found as (qq) / (qq) — 0.02 for Air — 0.5a~^ ~ 
0.4GeV around the physical region of m ~ 0.006a^^ ~ 5MeV. 



By removing the low-lying Dirac modes, the chiral 
condensate (qq) is largely reduced, reflecting the Banks- 
Casher relation. Actually, directly from lattice QCD 
calculation, we find a large reduction of the chiral con- 
densate in the presence of the IR cut Air = 0.5a~^ ~ 
0.4GeV, 



0.02, 



(36) 



around the physical region of m ~ 0.006a~^ ~ 5MeV 
as shown Fig. 2. 
Figure 3 shows the Dirac-mode projected inter-quark 
potential V^{R) after removing low- lying Dirac modes, 
which is obtained in lattice QCD with the IR-cut of 
Pir(A) = p(A)6'(|A| - Air) with the IR-cutoff Air = 
0.5a^^. Remarkably, no significant change is observed 
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1.5 



0.5 



where a certain part of Ai < |A„| < A2 of Dirac modes is 
removed. Unfortunately, when the wide region of Dirac 
modes is removed, the statistical error becomes quite 
large for the Dirac-mode projected inter-quark potential. 
Here, we remove the IM Dirac modes of 0.5 — 0.8[a^^], 
0.8 — 1.0[a~^], and 1.0 — 1.2[a~^], respectively, and inves- 
tigate the corresponding IM-cut inter-quark potential in 
each case, as shown in Fig. 5. For each case, the confining 
force is found to be almost unchanged by the IM-cut. 



R [a] 



FIG. 3: The lattice QCD result of the inter-quark poten- 
tial V^{R) after removing low- lying Dirac modes. The cir- 
cle symbol denotes the potential obtained with the IR-cut of 
Pir(A) = p(A)6l(|A| - Air) with the IR-cutoff Air = 0.5a"\ 
The square symbol denotes the original inter-quark potential. 
The potential is almost unchanged even after removing the 
low-lying Dirac modes, apart from an irrelevant constant. 



on the inter-quark potential besides an irrelevant con- 
stant, that is, the confining force is almost unchanged, 
even after removing the low-lying Dirac modes, which 
are responsible to chiral symmetry breaking [23|- This 
result indicates that one-to-one correspondence does not 
hold for confinement and chiral symmetry breaking in 
QCD. 

We also investigate the UV-cut of Dirac modes in 
lattice QCD, using the UV-cut Dirac spectral density 
puv(A) = p(A)6'(Auv - |A|) with the UV-cutoff Auv = 
2a~^ ~ 1.6GeV. In this case, unlike the IR cut, the chi- 
ral condensate is almost unchanged, and chiral symmetry 
breaking is almost kept. We show in Fig. 4 the UV- 
cut inter-quark potential, after removing the UV Dirac 
modes. We find that the confining force is almost un- 
changed by the UV-cut. This result seems consistent 
with the pioneering lattice study of Synatschke-Wipf- 
Langfeld [13 : they found that the confinement potential 
is almost reproduced only with low-lying Dirac modes, 
using the spectral sum of the Polyakov loop [s^l ■ 



FIG. 4: (a) The UV-cut Dirac spectral density puv(A) = 
p(A)6i(Auv - |A|) with the UV-cutoff Auv = 2a-^ ~ 1.6GeV. 
(b) The UV-cut inter-quark potential (circle) after removing 
the UV Dirac modes. The potential is almost unchanged from 
the original one (square), apart from an irrelevant constant. 



Furthermore, we examine "intermediate(IM)-cut" , 



I ^ 



1,5 2 25 



2 3 



X [a"'l 



1.5 2 2.5 




1.5 2 2.5 



FIG. 5: The left figures show the intermediate(IM)-cut Dirac 
spectral density pim(A): the IM Dirac modes of 0.5 — 0.8[a^^] 
(top), 0.8-1.0[a"^] (middle), and 1.0-1.2[a^^] (bottom) are 
cut. The right figures show the corresponding IM-cut inter- 
quark potential (circle) after removing the IM Dirac modes. 
For each case, the potential is almost unchanged from the 
original one (square), apart from an irrelevant constant. 



Thus, from the above lattice QCD results, we conclude 
that there is no specific region of the Dirac modes respon- 
sible to the confinement in QCD, unlike the chiral sym- 
metry breaking. Instead, we conjecture that the "seed" 
of confinement is distributed not only in low-lying Dirac 
modes but also in a wider region of the Dirac-mode space. 



VI. SUMMARY AND DISCUSSIONS 

We have developed a manifestly gauge-covariant ex- 
pansion using the eigen-mode of the QCD Dirac opera- 
tor p = 7''!?''. With this method, we have performed 
a direct investigation of correspondence between confine- 
ment and chiral symmetry breaking in lattice QCD. As 
a remarkable fact, the confinement force is almost un- 
changed even after removing the low-lying Dirac modes. 
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which are responsible to chiral symmetry breaking. This 
indicates that one-to-one correspondence does not hold 
for between confinement and chiral symmetry breaking 
in QCD. 

Our strategy is to investigate the relation between the 
nonperturbative properties of QCD, by extracting or re- 
moving the essence of chiral symmetry breaking. This is 
similar to the demonstration of Abelian/monopole dom- 
mance i^-ii, [11] or center/vortex dominance 

fisl l35l |36| for nonperturbative properties. However, 
while the previous scenario has been done in a specific 
gauge, our new method is manifestly gauge-invariant. 
In this analysis, we have carefully amputated only the 
"essence of chiral symmetry breaking" by cutting off the 
low-lying Dirac modes. Then, we have artificially realized 
the "confined but chiral restored situation" in QCD. 

Recently, Lang and Schrock studied the hadron spectra 
after the cut of the low- lying Dirac modes [s^ . Since the 
quark propagator is directly expressed with the Dirac op- 
erator the Dirac-mode projection is straightforward, 
and complicated projection procedure is not necessary 
in such studies. In their study, although the confine- 
ment was not checked, the appearance of hadronic spec- 
tra seems to suggest the existence of the confinement 
force, even after cutting the low-lying Dirac modes. 

Next, we comment on the possible relation among con- 
finement, chiral symmetry breaking, and monopoles in 
QCD. There is a close relation between confinement and 
chiral symmetry breaking through the monopoles in the 
MA gauge |12h14| . The monopolc would be essential de- 
grees of freedom for most nonperturbative QCD: confine- 
ment f2l|, chiral symmetry breaking and instan- 
tons [22|. In fact, removing the monopole would be "too 
fatal" for the nonperturbative properties, so that non- 
perturbative QCD phenomena are simultaneously lost by 
their cut. On the approximate coincidence of the critical 
temperatures of dcconfincmcnt and chiral restoration, a 
large change of monopoles may lead to both phase tran- 
sitions since the global connection of the monopolc 
current seems to be largely changed around the QCD 
phase transition j23j . 

As for the recent finite-temperature QCD analysis, a 
lattice QCD group has reported a certain difference be- 
tween the "critical temperatures" of dcconfinement and 
chiral restoration, which arc determined by the suscepti- 
bility peak of the Polyakov loop and chiral condensate. 



respectively [37[. This may be also an indirect evidence 
of "confinement ^ chiral symmetry breaking" in QCD. 

Next; we briefly discuss the role of low- lying Dirac 
modes in the viewpoint of instantons in QCD. The 
Dirac zero-mode associated with an instanton is local- 
ized around it [llj. However, the localized objects are 
hard to contribute to the large-distance phenomenon of 
confinement, although such low- lying Dirac modes con- 
tribute to chiral symmetry breaking. Recall that instan- 
tons contribute to chiral symmetry breaking, but do not 
directly lead to conflncment Then, as a thought 

experiment, if only instantons can be carefully removed 
from the QCD vacuum, confinement properties would be 
almost unchanged, but the chiral condensate is largely 
reduced, and accordingly some low-lying Dirac modes 
disappear. Thus, in this case, confinement is almost un- 
changed, in spite of the large reduction of low-lying Dirac 
modes. 

If the relation between confinement and chiral sym- 
metry breaking is not one-to-one in QCD, richer phase 
structure is expected in QCD. For example, the phase 
transition point can be different between dcconfinement 
and chiral restoration in the presence of strong electro- 
magnetic fields, because of their nontrivial effect on chiral 
symmetry [ssj . 

To obtain more definite conclusion, it is desired to per- 
form larger-volume lattice QCD calculations and to cut 
various wider region of the Dirac modes, although it is 
quite difficult. It is also interesting to investigate the sim- 
ilar analysis at finite temperatures in lattice QCD. The 
full QCD calculation in this direction is also an interest- 
ing subject. 
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